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ABSTRACT 

Context. Evolutionary trends in the surface abundances of heavier elements have recently been identified in the globular cluster 
NGC 6397 ([Fe/H]= -2), indicating the operation of atomic diffusion in these stars. Such trends constitute important constraints for 
the extent to which diffusion modifies the internal structure and surface abundances of solar- type, metal-poor stars. 
I i Aims. We perform an independent check of the reality and size of abundance variations within this metal-poor globular cluster. 

Methods. Observational data covering a large stellar sample, located between the cluster turn-off point and the base of the red 
giant branch, are homogeneously analysed. The spectroscopic data were obtained with the medium-high resolution spectrograph 
FLAMES/GIRAFFE on VLT-UT2 (R ~ 27 000). We derive independent effective-temperature scales from profile fitting of Balmer 
lines and by applying colour- T e g calibrations to Stromgren uvby and broad-band BVI photometry. An automated spectral analysis 
code is used together with a grid of MARCS model atmospheres to derive stellar surface abundances of Mg, Ca, Ti, and Fe. 
Results. We identify systematically higher iron abundances for more evolved stars. The turn-off point stars are found to have 0. 13 dex 
lower surface abundances of iron compared to the coolest, most evolved stars in our sample. There is a strong indication of a similar 
trend in magnesium, whereas calcium and titanium abundances are more homogeneous. Within reasonable error limits, the obtained 
abundance trends are in agreement with the predictions of stellar structure models including diffusive processes (sedimentation, 
levitation), if additional turbulent mixing below the outer convection zone is included. 

Key words. Stars: abundances - Stars: Population II - Globular clusters: general - Techniques: spectroscopic - Methods: observa- 
' tional - Diffusion 
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C^) 1. Introduction High ages and thin convective envelopes are factors that make 

I " ■ metal-poor TO stars particularly prone to atomic diffusion. 

> A globular cluster (GC) constitutes, in many respects, a homoge- Qnce a star leayes the MS and begins tQ deyelop Mq ±e 

^ ; neous stellar population and is therefore an excellent laboratory iant ks QUter convection zone gradua n y reac hes deeper 

' f ° r teStmg and constrainm g models of stellar stmcture and layers, and elements that had previously been drained from 
d ■ evolution. Chemical abundances inferred from observations of ±e surface layers are mixed up again> eventually restoring the 

GCs can give important clues to the physical processes at work chemical composition (with fragi i e elements such as Li as 

in individual stars. Recently a discovery was made pointing to nQtable exceptions) . This is predicted to result in a gradual rise 

the existence of systematic differences in surface abundances in ±e surface abundances of the depleted elements as the star 

between stars m different evolutionary stages in a metal-poor evolyes al ±e sub iant and red iant branches (her e a fter, we 

GC. | Kornetal. | <|2007j) (hereafter Paper I) observe four samples write SGB for subgiant branch) GC ^ presumably shal - e the 

of stars located between the main sequence (MS) turn-off (TO) same original composition of elements such as iron, titanium, 

point and the red giant branch (RGB) of NGC 6397, using and which fa combination with their hi h and low 

FLAMES/UVES on the VLT, and found significant variations metal content make them suitable tefjt cases for stellar stmctme 

in abundances with effective temperature and surface gravity. modelfj including atomic diffusion . 

In particular, they conclude iron to be under-abundant by The abundance trend fo r koQ presente d in Paper I is directly 

0.16 ± 0.05 dex in the atmospheres of the TO stars compared to contrad i ctory t0 result bv lGratton et all (flOOlh . who by studying 

the RGB stars in their sample. Proposed as an explanation to this TQ ys base . RGB (bRGB) stars conc i ude that there are no 

phenomenon is a continuous, long-term depletion of iron and si ificant diff e r ence in iron abundance between the groups 

other heavy elements from the surface layers, at work during (magnesium ^ howeV er, found to be 0. 15 dex less abundant in 

the stars' life-time on the MS. The mechanism responsible the TO stars). The cause for the differing result is twofold: lower 

for the depletion is generally referred to as atomic diffusion - effectiye te matures are found for the TO stars in Paper I and 

an umbrella term accounting for several diffusive processes. ±e analysifj indude starfj Qn ±e RQB 

~ Here we present a follow-up analysis of spectroscopic data 

Based on data collected at European Southern Observatory (ESO), , pr AMFS/GIRAFFF collected simultaneouslv to the 

Paranal, Chile, under program ID 075.D-0125(A). tr0m ^AMfc,J>/UiKAr<r<fc, collected simultaneously to tne 
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observations presented in Paper I. A large sample of stars, 
~ 100, covering the range from the TO to the bRGB is ho- 
mogeneously analysed, with the aim of further constraining 
possible variations in surface abundance in the cluster. In 
particular, we investigate whether the results presented in Paper 
I are robust for a larger sample of stars, analysed with an 
independent set of tools. We start by describing the observations 
and the data-reduction procedure in Sect. 2. Sect. 3 is devoted 
to the determination of fundamental stellar parameters and a 
description of the spectrum analysis code. In Sect. 4 we present 
and discuss our abundance results. Sect. 5 compares predictions 
from stellar-structure models including atomic diffusion to the 
abundances we have inferred from observations. In Sect. 6 we 
summarise our conclusions. 

2. Observations and data reduction 

2.1. Photometric observations 

The target selection for the spectroscopic study is based on 
Stromgren uvby photometry. The photometric observations 
were collected with the DFOSC instrument on the 1.5 m tele- 
scope on La Silla, Chile, in 1997. Additional BVI photometric 
data were obtained in 2005. The UCAC2 catalog at Vizier 
dOchsenbein et al.l l2000i) provides astrometry with a precision 
of ~ O'.'IO for individual stars. A full description of the photo- 
metric observations and their reduction is given in Paper I and 
references therein. 

We select a sample of 135 stars for spectroscopic analysis, 
located between the cluster TO point at V — 16.3 to the bRGB 
at V = 15.2. Apart from providing the target list, the photometry 
is also considered in this work in the determination of stellar 
parameters. Fig. Q] shows a colour-magnitude diagram, V versus 
(v - y), of the cluster, with the selected stars marked. 



2.2. Spectroscopic observations 

All spectroscopic data were collected in Service Mode, with 
the fibre-fed, multi-object, medium-high resol ution spectro- 
graph FLAMES/GIRAFFE dPasquini et al.ll2002l) at ESO-VLT. 
FLAMES allows for 132 objects to be observed simultaneously, 
with GIRAFFE in MEDUSA mode. A total exposure time of 
20.5 h was spent under grey-to-bright lunar conditions, but with 
good seeing (average seeing ~ 0.7"). We apply four of the high- 
est resolution 'B-settings', which together cover a total spectral 
range of 822 A in the optical. On average 15 fibres per observa- 
tion are dedicated to a simultaneous monitoring of the sky back- 
ground to ensure a proper sky-correction. Stars in different evo- 
lutionary stages are distributed randomly across the FLAMES 
field-of-view. We find no gradient in background light over the 
field and subtract a single, averaged sky-spectrum from the stel- 
lar spectra obtained during each observation. The setup is briefly 
summarised in Table Q] 

Basic reduction of the spectroscopic data is performed with the 
ESO-maintained GIRAFFE pipeline, version 1.0, using stan- 
dard settings. Further processing of the data is carried out with 
MIDAS and C-routines (kindly provided by N. Christlieb). The 
GIRAFFE CCD chip has a defect in its upper right corner caused 
by the illumination of a diode, giving elevated dark-current val- 
ues to a significant number of pixels. If uncorrected for, this 
will propagate into a flux upturn towards longer wavelengths 
in dozens of spectra. The feature can be removed by subtract- 
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Fig. 1. Observed colour-magnitude diagram of NGC 6397. The 
targets selected for the spectroscopic study are marked with 
black squares. 



Table 1. Observational setup for FLAMES/GIRAFFE ". 



Setting 


Wavelength 


Resolution 


Exposure 


Dates 2005 


# 


range [A] 


A/AA 


time 




5B 


4376-4552 


26000 


7h40m 


23,26-29 Mar 


7B 


4742-4932 


26 700 


5hl0m 


29-31 Mar 


9B 


5 143 - 5 356 


25 900 


3h50m 


01-03 Apr 


14B 


6 383 -6 626 


28 800 


3h50m 


02-04 Apr 



" Data from http://www.eso.org/instruments/flames 



ing a dark frame from the raw science frames (after proper 
scaling with exposure time). To avoid the introduction of addi- 
tional noise we chose to correct for dark current after the spectra 
had been extracted, when the glow is a well-behaved and easily 
smoothable function of wavelength. For this purpose we pseudo- 
reduce a dark frame produced close in time to the observations. 
We find, however, that it is necessary to apply an additional, mi- 
nor scaling factor to the dark counts to completely remove the 
feature, which we ascribe to the documented time-variability of 
the glow strength^ 

The radial velocity correction of our stellar spectra results in 
a cluster mean heliocentric radial velocity of 18.1 ± 0.3kms _1 
from measurements of 133 stars. We find a velocity dispersion 
of 3.40 kms -1 . This radial velocit y value is in agreement with 
the recent result by iMilone et al.l (120061) who derive 18.36 + 
0.09(+0.10)kms _1 (the first is the random error and the second 
the systematic error), obtained from 1486 stars. We identify two 
stars that are not members of the cluster, 502219 and 501856, 
both with deviating radial velocities and evidently too metal-rich 
spectra (see Sec. 4). 



1 see the instrument's quality control pages: 
http://www.eso.org/observing/dfo/quality/GIRAFFE/qc/dark_qcl.html 
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Fig. 2. The best-fit Ha profiles for three stars. The thick line cor- 
responds to a model with the stated effective temperature. The 
grey-shaded regions represent the regions used for fitting. 

After sky-subtraction and radial velocity correction the spectra 
are coadded and rebinned onto a final wavelength scale with a 
step size of 0. 1 A. The S /Af -ratio of the coadded spectra in the 
14B setting varies from 28 to 148 per rebinned pixel with a mean 
value of 85, estimated from a relatively line-free region between 
6 440 -6 449 A. The corresponding mean value for the bluer 7B 
setting is 76, from measurements between 4812-4 821 A. 

3. Analysis 

In this section we describe how we derive the fundamental pa- 
rameters necessary to conduct the abundance analysis. We em- 
phasise that we have made an effort to carry out all steps of 
the investigation using independent tools from those of Paper 
I, in particular by implementing different spectral synthesis and 
model atmospheres codes. 

3.1. Effective temperatures 

Profile fitting of Balmer lines is a well-established r e ff-indicator 
in solar-type stars. Considering the reliable flat-fielding which 
can be obtained with fibre-fed spectroscopy and the large free 
spectral range of GIRAFFE (FSR ~ 200 A), this data set is well 
suited for analysis of broad lines such as the Balmer lines. We 
implement an automated ^-minimization technique for the nor- 
malisation and profile fitting of the wings of Ha. Synthetic 
spectra are produced with the FORTRAN code HLINPROF 
(Barklem & Piskunov 20030, which is incorporated in the LTE 
(local the rmodynamical equilibrium) spectrum synthesis code 
SYNTH dPiskunovll 1992ft . We t reat Stark broadening acco rding 
to the tabulated fine profiles of IStehle & Hutcheonl d 1999ft and 
self broadening (res onance broadening an d van der Waals broad- 
ening) according to Bar klem et al.l d2000l) . Radiative broadening 
and an estimate of line broadening due to helium collisions are 
also included in the calculations. 

To compute a set of plane-parallel model atm ospheres for the 
spectrum synthesis we use the MARCS code (Gustafsson et alj 
ll975MAsplund et al.lfl9"97ft . MARCS is a one-dimensional (ID), 
LTE model atmospheres code, invoking flux conservation (ra- 

2 Available at www.astro.uu.se/~barklem. 
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Fig. 3. Comparison between our Ha-based effective tempera- 
tures and those obtained by applying the calibrated relations of 
lAlonso et al. I (1 19961 1 19991) to colour indices v-y and V-I. The 
y-axis shows Ar e ff = r e ff,Ha - ^photometry for the stated index and 
calibration. 

diative and convective) and hydrostatic equilibrium. Convection 
is treated according to the customary mixing length theory, with 
the mixing length parameters set to a - 1 .5 and y = 0.076. 
The observed spectrum is compared to a grid of synthetic spec- 
tra, equidistant in effective temperature with 10 K steps. The fit- 
ting only considers regions that are not significantly affected by 
absorption from metal and telluric lines (see Fig. |2]l as judged 
from high-resolution and high-S /N spectra of metal-poor sub- 
giants. First, an initial guess of the correct model is made and 
only wavelength regions located +35 - 50A from the Ha-line 
center are used. Here the sensitivity of the hydrogen line to ef- 
fective temperature is low and the difference between observed 
and modelled flux is minimised to find the optimal straight line 
used for normalisation. Thereafter, only wavelength regions in 
the effective-temperature sensitive parts are used, located within 
+25A of the center, to find the model that minimises the x 2 ~ 
distance to the observed flux. However, to avoid modelling un- 
certainties affecting the line core, only residual flux above 0.75 
is considered. The whole procedure is iterated until the initial 
guess coincides with the resulting r er f-value. 
By repeatedly adding random noise to synthetic spectra and 
feeding them as input to the procedure, the precision of the fit- 
ting method can be determined. We find that a typical S /W-ratio 
of 85 gives a random error of about 30 K, for the coolest as well 
as for the hottest stars. For our lowest-quality spectra the corre- 
sponding value is 80 K and for the highest quality data the ran- 
dom error is as low as 15 K. The absolute error in T e ff is more 
difficult to assess, as it depends on several factors in addition 
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to observational uncertainty, of which the most important ones 
are the line broadening of hydrogen, the treatment of convection 
in the s tellar atmosph e re, an d the stellar parameters logg and 
[Fe/Hl . iBarklem et al.l (|2002) show that these contributions to- 
gether may add up to an error of the order of 100K for stars in 
this metallicity and effective-temperature range. 
From the FLAMES solar atla^3 we retrieve reduced GIRAFFE 
spectra from observations illuminating all MEDUSA fibers in 
the 14B-setting, which we average into a single solar spectrum. 
Using the setup described above we derive T^^a = 5 630 K for 
the Sun, which is approximately 150 K lower than the accepted 
value, an offset which we as cribe mainly to observational uncer- 
tainty. Barklerr Tet al.l (120021) derive a solar effective temperature 
of 5 722 + 8 1 K with th e same models, us ing an Ha-spectrum 
from the Kitt Peak atlas (Ku rucz et al.ll 984). The Ha line profile 
found from the FLAMES/GIRAFFE solar atlas is more narrow 
compared to e.g. the Kitt peak atlas and we recommend to re- 
observe the Sun with GIRAFFE with the new CCD (upgraded in 
May 2008) and seek the reason for this difference. 
Using our obtained solar effective temperature as a zero-point 
offset, we consequently shift all effective temperatures by 
+ 150K. We thus assume that, in the stellar-parameter range 
spanned by our targets, the Ha-line has a similar sensitivity to 
effective temperature as the Sun, which may be realistic to a first 
order approximation. The shifted Ha-based effective tempera- 
tures are the ones we later adopt in the abundance analysis and 
we will refer to them as the spectroscopic T e ff in the remaining 
part of the paper. 

To assess the validity of our spectroscopic effective temperatures 
we construct photometric r e ff scales fr om four different colour 
indices. The lAlonso et al.l (1 19961 11999) relations, calibrated on 
the infrared flux method, are used to derive effective tempera- 
tures based on b - y, v - y, B - V, and V - I. Note that, at this 
metallicit y and stellar-parameter spac e, the corresponding rela- 
tions of Ramr rez~& Melendezl d2005l) produce very similar ef- 
fective temperatures (see Paper I). There are two Alonso et al. 
calibrations for each inde x, one suitable for main sequence stars 
of spectral types F0-K5 dAlonso et al. 1996) and o ne for giant 
stars in the same spectral range dAlonso et alJll99 9Fl Since our 
targets range from TO stars at the very end of the MS to bRGB 
stars, both the dwarf and giant calibrations are used to separately 
derive effective temperatures. Note that by doing so, we also ap- 
ply the T e ff - colour relations to stars that do not fall in the range 
in colour covered by the calibrations at this metallicity, which 
may affect the results. 

Before applying the calibrations, we eliminate star-to-star scat- 
ter caused by errors in the observed colours. We clean the 
photometric sample based on data quality, as measured by the 
DAPHOT SHARP parameter, and construct fiducial relations 
between each colour index and the V magnitude, which is the 
best observed and calibrated quantity. The relations are obtained 
by averaging the colours for stars in V bins of 0.22 magni- 
tudes. We then shift the observed colour for each star onto the 
sequence and, finally, we correct all colours for the reddening 
of NGC 6397. This value has been measured in several studies 
and typic al estimates l ie in the range E(B - V) = 0.17 - 0.20 
(see e.g. Reid & Gi zis| 1998). Here a value o f 0.179 is adopted, 
following Anthonv-Twarog & Twarog (2000). The reddening in 
the other colour indices are derived from E(B - V), using the 



3 http://www.eso.org/observing/dfo/quality/GIRAFFE/pipeline/ 
solar.html 

4 Calibrations for the v - >'-index for dwarf and giant stars are pub- 
lished in Paper I. 



relation coefficients given in Ramirez & Melendez ( 2005]) and 
E(y-y) = 1.7 xE(b-y). 

Fig. [3] shows a comparison between our derived spectroscopic 
and photometric effective temperatures for 122 stars (all stars for 
which the Ha line is observed). The difference between the Ha- 
based values and the effective temperatures obtained from the 
dwarf and giant calibrations of the narrow-band index v — y and 
the broad-band index V - 1 are plotted against T e sp a . The agree- 
ment between Ha and the v—y dwarf calibration scale is good for 
the hotter half of the sample, with spectroscopic T e s > 5 900 K. 
The average difference for individual stars is here ~ 70 K. For 
the cooler half, at spectroscopic T e s ^ 5 900 K, good agreement 
is instead seen with the giant calibration, with a mean difference 
~ 50 K. The spectroscopic temperature scale hence suggests a 
transition between the two v — y calibration scales at some point 
in the middle of the SGB, in line with what may be expected. 
The other narrow-band index, b - y (not shown), produces very 
similar effective temperatures to the v—y scales. The V — I dwarf 
calibration scale appears similar in relative behaviour to v - y, 
but the effective temperatures are shifted towards slightly cooler 
values, giving a good agreement with the spectroscopic T e ff for 
the cooler half of the T e ff-range. The giant scale, however, is ev- 
idently too cool overall to reproduce the spectroscopic results. 
For B - V (not shown) the two Alonso et al. calibrations result 
in scales that are identical to each other to within 5 K. In com- 
parison to the Ha-based values the agreement is reasonable at 
the hotter end whereas the cooler end is offset to higher tem- 
peratures by up to ~ 200 K. This index thus implies a markedly 
shorter total effective temperature range for our sample than the 
other three. 

As an additional test, we also compared the effective tem- 
peratures obtained with the Alonso et al. relations to (v - y) 
and (b - y)-ca librations based on synthetic MARCS colours 
dGustafsson et alj 2008; calibrations to be published in Onehag 
et al. 2008, in preparation). The theoretical calibrations repro- 
duce the r e ff-range spanned by the targets well, indicating a 
somewhat shorter total range, by about ~ 70 K, than what is pre- 
dicted from a combination of the two Alonso et al. relations for 
b — y and v — y, respectively. 

Much effort has previously been made to pin down the effective- 
temperature scale of this cluster at the SGB. A summary of the 
results from six high-resolution studies, including Paper I, is 
given in Fig. [4] along with our spectroscopic and four photo- 
metric scales. In the plot, which illustrates visual magnitude ver- 
sus effective temperature, filled symbols represent previous re- 
sults, and lines represent this study. For clarity, we show only a 
smoothed relation between V and T^Ha (the lcr-scatter around 
this line is 90 K) and only parts of the phot ometric V - r e ff re- 
lations based on v - y and V - I, i.e. the Alo nso etal.1 dl 996) 
calibra tion is applied only to the hotter half and the lAlonso et all 
d 19991) only to the cooler half of the stellar sample. All data are 
collected from the stated publications or from references therein. 
Only the effective temperatures that are actually adopted by the 
different authors are shown, even if several different scales are 
presented in the study . The visual magnitudes in the sample of 
iThevenin et a fl d2001l) are from their Table 1 , taking the average 
whenever two values are given. 

Fig. [4] shows that the different studies agree within approxi- 
mately 200 - 300 K at a given visual magnitude, which may 
be a reasonable error, albeit undesirably large when the aim is 
precise chemical abundances. Most studies focus on TO stars 
and the only two studies, except the pres ent, that include star s 
on both ends of the SGB are Paper I and lGratton et alj d2001l) . 
On a relative scale, which is most essential for this work, we 
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Fig. 4. Comparison between our spectroscopic r e ff-scale, the 
photometric values obtained from v-y and V — I, and the re- 
sults from six other studies. 



predict a slightly larger effective temperature range between 
V = 15.5 - 16.2 than Paper I and a shorter range than Gratton 
et al. In the former study, effective temperatures are also ob- 
tained by Ho- profile fitting but using a different model atmo- 
spher e and spectrum synthesis code (ODF-MAFAGS, Grupp 
2004, Fuhrma nn et al.l 1993). For the eight targets in common 
with this study we find group-averaged effective temperatures 
that are different by less than 10K. The typical difference for 
individual stars is 50 K. The T e ff value found for the TO-star 
group of the Gratton et al. study exceeds our estimate by almost 
200 K and is extensively discussed in Paper I, where the authors 
point out that the effective-temperature assignment of Gratton et 
al. , based on Ha, may be systematically affected by the echelle- 
order blaze residuals imprinted on the continuum in the UVES 
(slit-mode) observations. Differences in model atmospheres and 
hydrogen line-broadening theory used may also explain part of 
the offset (see Paper I). 

3.2. Surface gravities 

We derive surface gravities for the target stars using the custom- 
ary relation between effective temperature, luminosity, mass, and 
surface gravity. Luminositie s are calcula t ed fro m the apparent 
visual magnitude V and the lAlonso et al.l {T999) calibration for 
bolometric correction, which is given as a function of metallic - 
ity, [Fe/Hfl and r e ff. The metallicity is set to -2.0 for all stars. 
The distance modulus of NGC 6397 is assumed to be 12.57. 
Ste llar masses are inferr ed from a 13.5 Gyr isochrone of the clus- 
ter dRichard et al.l2005l) . which places the stars in the mass range 
0.78 Mq -0.79 M0. With these values, the surface gravity found 
for TO stars at r eff = 6 250 K is log g = 3 .96 and log g = 3 .40 for 
bRGB stars at T e s = 5 450 K. We use a single, averaged value of 
log g for all stars for which we have derived the same T e g . 
The aim of this study is to draw conclusions about abundance 
differences between stars. We are therefore mainly interested in 
the accuracy with which we can determine surface gravities on 
the relative, rather than the absolute scale. In this respect, the 
effective temperatures have the largest, albeit small, influence 
on the surface gravity values. A rise in effective temperature of 



Fig. 5. A selected part of two observed, normalised spectra (dot- 
ted lines) and the best fitting model spectra (solid lines), as deter- 
mined by the analysis code. Note that the models in this figure 
only contain the spectral lines that are used in the abundance 
analysis. These lines are marked with labels. The top spectrum 
belongs to the bRGB star 23267 and the bottom plot to the TO 
star 12318 



100 K corresponds to an increase in logarithmic surface grav- 
ity of approximately 0.03 dex. This can be compared to an in- 
crease in stellar mass by 0.01 M0 that propagates into a rise in 
log g by 0.005 dex. We set a constant metallicity, to avoid cir- 
cular arguments as regards the existence of abundance trends. 
The influence of [Fe/H] on the derived surface gravities is any- 
way small. An increase in [Fe/H] by 0.1 dex propagates into a 
decrease in logg of 0.006 dex. We estimate the relative error in 
V to approximately 0.01 magnitudes, which gives a contribution 
to the error in \ogg of 0.004 dex. The value of the distance mod- 
ulus only affects the absolute values of the surface gravities. An 
increase in (m - M)v of 0.5 corresponds to an overall rise in 
logg of 0.2 dex. Given these low dependencies of relative pho- 
tometric surface-gravity values on other parameters, abundance 
differences within the sample can be constrained from gravity- 
sensitive lines (e.g. those of Fell) with high confidence. 



3.3. Spectrum analysis code 

For abundance analysis of the stellar sample we use the same 
automated code as developed for t he Hamburg/E SO R-process 
enhanced star survey, (HERES, Barklemetal] 2005). We 
summarise here som e important asp ects of the analysis code, 
referring the reader to Bark lem et al. 



for a full description. 



5 [X/Y] 
element X. 



log - log (jfy)o where N\ is the number density of 



The software is built on the I DL and C++ based Spectro scopy 
Made Easy (SME) package by IValenti & Piskunovl (U996). with 
subsequent updates, but has been modified to run without any 
user interaction. Some improvements have also been made to 
the original SME spectrum synthesis code, most importantly the 
inclusion of continuous scattering in the source function. As the 
software is implemented, we supply it with 1 ) a count spectrum 
and the measurement error in each pixel, 2) the final stellar 
parameters T e ff and log g and an initial guess for the metallicity 
([Fe/H]= -2.0), 3) a list of lines to base the abundance analysis 
on and wavelength w indows where each l ine can be found and 
4) a grid of MARCS (lAsplund et al.ll 19971) plane-parallel model 
atmospheres, all with scaled solar abundances except for an 
enhancement of alpha-elements of 0.4 dex. This setup is the 



6 



K. Lind et al.: Atomic diffusion and mixing in old stars 



same as implemented in Barklem et al.l (120051) . except that we 
keep log g fixed. 

The code automatically normalises a spectral region stretching 
7 A on both sides of each given wavelength window, by iterative 
fitting of a low-order polynomial. At this point any remaining 
cosmic ray hits are also identified and discarded. The line central 
wavelength is determined and, if necessary, shifted slightly 
(within the estimated uncertainty) to coincide with the model 
wavelength. Observed and modelled spectra are then compared 
and the best statistical match is de t ermin e d via a parameter - 
optimization algorithm (Marquardt] 119631: iPress et ail Il992l) . 
We assume the stars to be slow rotators and fix the projected 
rotational velocity to lkms -1 . The free parameters are thus 
the abundance, the microturbulence and the macroturbulence 
Vmacro (Gaussian, isotropic), in which also the instrumental 
broadening is included. The online appendix contains a list of 
all lines used in the abundance analysis, with references to the 
^/-values that were adopted. Based on the availability of lines 
covered in the spectra, we choose to derive abundances of iron, 
titanium, calcium, and magnesium. Fig. [5] shows a selected 
part of two observed, normalised spectra and the best-fit model 
spectra. 

The analysis assumes LTE line formation in one-dimensional, 
plane-parallel, model atmospheres. These simplifying assump- 
tions can, of course, influence the derived abundances to a 
certain extent, which we discuss further in Sect. 4. In particular, 
strong line cores are expected to be poorly described by LTE 
photospheric models. To avoid large biases, pixel values below 
50 % of the continuum flux are generally disregarded. 
To analyse the resulting abundances, we also need knowledge 
about how errors in the stellar parameters affect them. By rerun- 
ning the code with one perturbed parameter at a time, holding 
the others fixed, the abundance sensitivity to these quantities 
are found for the elements considered. We make three such ad- 
ditional runs, with r e ff + 100K, log g +0.3 dex, and £+0.2 km s _1 . 
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Fig. 6. A histogram of the slopes obtained for linear fits of iron 
abundance with excitation potential of Fe I lines. 



5 450 K and 6 250 K, and the estimated lcr scatter in log The 
microturbulence values are found to gradually decrease from 
a value 1.86 + 0.13 km s~' for the stars at = 6 250K to 
1 .47 + 0.09 km s -1 for the stars at 5 450 K. 
A preliminary analysis of the stars found to be non-members 
based on radial-velocity measurements (see Sect. 2) reveals that 
they are both subgiants, with approximate stellar parameters 
r eff = 5 700K and logg = 3.3. We find [Fe/H] = -0.6 
for 501856 and [Fe/H] = -1.1 for 502219, both indicating 
[Ca/Fe] = 0.3. Among the cluster stars, 7720 is found to 
have consistently outlying abundances. With T e s = 5 580K, it 
displays abundances exceeding the average by 0.19 - 0.26 dex. 
The star has excellent agreement between Ha-based and 
photometry-based effective temperatures and a radial velocity 
typical of the cluster. 



4. Results 

Fig. [7] shows our derived chemical abundances of magnesium, 
calcium, and iron, for 116 stars. Six stars, for which we have 
only two observations each, were rejected due to largely 
deviating abundance results (not shown), most likely caused by 
the low quality of these spectra. General trends are obtained by 
averaging the abundances in effective-temperature bins of 500 K 
(with shrinking box-size towards each end) and then applying 
additional smoothing to avoid the influence from outliers. The 
thin lines shown in the plot indicate the standard deviations of 
the trends' residuals. For comparison, Fig. [7] also shows the 
results from Paper I (open squares). Note that the abundances 
of Fe, Mg, and Ca then were derived by non-LTE (hereafter 
NLTE) analysis. 

It is evident that our derived abundances indeed suggest 
that there are trends with effective temperature on the SGB. 
Especially in iron there is a notable trend of increased abun- 
dances towards lower T e g. The magnesium values display 
a similar behaviour, but due to larger star-to-star scatter the 
trend is less clearly visible. The calcium abundances may 
point to a slight, insignificant, increase towards cooler ef- 
fective temperature, whereas the titanium abundances rather 
indicate the opposite, equally insignificant, trend with T e ff. 
Table [2] summarises the average abundances we derive for all 
elements at two representative effective temperature points, 



4.1. Error sources 

There are several possible sources of error to consider, which 
may effect our abundances on the relative and absolute scale. 
First of all our analysis relies on the assumption of LTE line 
formation, which is not in general a good approximation 
for metal-poor stars. Particularly lines from neutral minority 
species are believed to form out of LTE. NLTE corrections 
are expected to be larger in cool metal-poor stars than in their 
metal-rich counterparts, mainly because the reduced amount of 
metals gives less line-blocking in the UV, which causes more 
over-ionisation. Also, the amount of free electrons available in 
the atmosphere is reduced in metal-poor stars, which means less 
LTE-establishing electron collisions. Considering the narrow 
range in stellar parameter space covered by our targets, we 
would expect NLTE corrections to be similar for all stars, 
thereby mainly influencing the abundances on the absolute 
scale. 

The second point we consider is the abundance dependence on 
stellar parameters, most importantly on effective temperature. 
Table|3]shows what effect perturbing T e ff, log g, and £ from their 
derived values has on the derived abundances of all elements, 
at two effective temperature points. With this knowledge it is 
easily predicted how systematic errors in e.g. our r e (j scale 

6 log* = log(&) + 12 
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Fig. 7. Abundance trends (filled dots) with effective temperature for Mg, Ca, Ti, and Fe from analysis of 116 stars. The thick lines 
display general trends and the thin lines show estimated Icr limits. Open squares represent the results from Paper I, in which Mg, 
Ca, and Fe are treated in NLTE and Ti in LTE. 

Table 2. Average abundances obtained at two effective temperature points. 



7-eff [K] 


logg [cgs] 


f [kms H ] 


log £Mg 


0"Mg 


log 6 Ca 


<7"Ca 


log en 


<7"xi 


log 6 Fc 




6 250 


3.96 


1.86 


5.37 


0.10 


4.31 


0.06 


3.03 


0.05 


5.09 


0.07 


5 450 


3.40 


1.47 


5.54 


0.10 


4.38 


0.06 


3.00 


0.05 


5.22 


0.07 


A 800 


0.56 


0.39 


-0.17 




-0.07 




0.03 




-0.13 





would influence the results. When discussing this issue for 
each element, we take into account that changes in the effective 
temperature scale cause changes in surface gravity (see Sect. 
3.2). We neglect possible changes in £ caused by altering T e g or 

The most influential parameter on the chemical abundances 
is the effective temperature. As seen in Sec. 3.1, T e ff- 
determinations based on Ha and cluster photometry seem 
to be in good agreement, but additional information can be 
drawn from the excitation equilibria of the targets. Generally, 
the abundance of each element is the one that gives the optimal 
match to all lines, but for testing purposes we also compute 
individual abundances for all neutral iron lines and thus can 
compare our adopted temperature scale with that implied by 
excitation equilibrium. The resulting slopes obtained for iron 
abundance with excitation potential of the lines are shown in a 
histogram in Figure [6] The median slope for the sample is neg- 
ative, -0.038 dex/eV, i.e. overall cooler effective temperatures 
are needed to establish the equilibria. By estimate, the equilibria 
of the hottest stars, T e g > 6000 K, are reproduced with a 



lowering of the effective temperatures of approximately 270 K. 
The stars in the cooler half, T e g < 5800 K, require a similar 
change, 260 K, as the abundance slopes with excitation potential 
are more negative than the hotter half, but also have a higher 
temperature-sensitivity. There is thus no sign of a sizeable 
expansion or contraction of the adopted effective-temperature 
scale. 

A third source of error originates in the spectrum itself. We 
may judge how well GIRAFFE performs in comparison to the 
higher resolution spectrograph UVES as eight of our observed 
stars, five TO stars and three bRGB stars, are the same targets 
as presented in Paper I. The results from measurements of a 
number of equivalent widths of lines in the GIRAFFE and 
UVES spectra of these stars are displayed in Fig. [8] The 
equivalent widths measured in the GIRAFFE spectra are 
plotted against the corresponding value obtained from the same 
lines in the UVES spectra. Filled circles represent TO stars 
and open circles bRGB stars. The results indicate that the 
GIRAFFE spectra have slightly weaker lines compared to those 
obtained with UVES. The difference amounts to approximately 
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Fig. 8. The figure shows a comparison between equivalent 
widths measured for the same lines in spectra of the same stars, 
obtained with GIRAFFE and UVES, respectively. The filled cir- 
cles corresponds to stars located at the TO and the open circles to 
bRGB stars. The dashed line mark 1:1 agreement and the solid 
line is a linear fit to the measurements. 



4%, which in logarithmic abundance units corresponds to 
-0.02 dex (assuming a linear relation between abundance and 
equivalent width). There is no obvious difference between 
TO and bRGB stars, leading us to conclude that this has no 
impact on the relative abundances we derive for our sample stars. 



4.2. Magnesium 

Our magnesium abundances are based on two of the Mglb 
triplet lines, 5 172 A and 5 183 A. As seen in Table [3] mag- 
nesium is the most sensitive of all elements to T e s and logg, 
which may contribute to explaining the higher level of scatter, 
approximately 0.10 dex, seen in logfMg compared to the other 
elements. Another probable source of scatter is conversion of 
Mg to Al, which is observed as an i nverse correlat ion between 
the two elements in GC stars (see e.g. lGratton et al.l 2004). 
As seen in Table [2] we obtain a difference of 0.17 dex between 
stars at the cooler and hotter end of the SGB. For a trend of 
this size to vanish, i.e. for our magnesium abundances to agree 
over the whole sample, the effective temperature of the TO 
stars would have to be raised by at least 320 K, and more, 
around 450 K, if one takes into consideration a simultaneous 
rise in surface gravity of 0.03 dex per 100 K (See Sect. 3.2), 
since the two stellar parameters have opposite influence on the 
magnesium abundance. Alternatively, the effective temperature 
of the bRGB stars would have to be lowered by 170 - 220 K. 
We note that the trend we obtain of 0.17 dex is larger than 
the 0.08 dex found in Paper I for the difference in magnesium 
abundance between their bRGB stars at average r e jf = 5 456 K 
and TO stars at r eff = 6 254 K. 

Based on existing studies, NLTE corrections for lines of 
neutral magnesium are positi ve in these st a rs, du e to over- 
ionization. Paper I, following iGehren et al.l (120041) . calculate 
logeNLTE - log 6lte to range between 0.11 dex and 0.1 4 dex 
in their sample, based on Mg I 5 528 A. Hence there are only 
small differential e ffects in abundance, judging from this line. 
Mas honkinaetalJ d2008l) find an average NLTE correction of 



0.12 dex for the somewhat hotter, but otherwise very similar, 
metal-poor star HD 84937, from five Mgl lines, including 
5 528 A, 5172 A, and 5 183 A. 

NLTE effects are pronounced in the line cores of the triplet 
lines, in this case leading to the observed lines being deeper 
than the LTE-based prediction. These influences are partly 
avoided for the coolest stars in the sample, for which the 
line cores extend deeper than 50% of the continuum flux and 
are disregarded by the spectrum analysis code. In the hottest 
stars, all pixels contributing to the lines are evaluated. Given 
this methodology, differential effects on the abundance trend 
obtained from photospheric NLTE modelling, including the line 
core, cannot be excluded. 



4.3. Calcium 

Nine lines of neutral calcium are considered in the determi- 
nation of log eca- The strongest lines, i.e. those consisting of 
the highest number of pixels and thereby having the largest 
influence on the derived abundances, are 4 434 A, 4 454 A, 
6 462 A, and 6 439 A. Our average calcium abundances increase 
from 4.31 dex at the TO to 4.38 dex at the bRGB, i.e. a trend 
of 0.07 dex. This can be compared to 0.02 dex found in Paper 
I, based on 6 122 A, 6 162 A, and 6 439 A. Table [3] shows that 
log eca has a low sensitivity to changes in T e ff and log g and the 
scatter in abundance is also small. To completely remove a trend 
of 0.07 dex in calcium abundance, one would have to raise the 
TO star effective temperatures by 500 - 750 K (assuming the 
dependence is linear also for such a large perturbation), again 
depending on the corresponding change in surface gravity, or 
equivalently lower the bRGB star effective temperatures by 
210- 260 K. 

Following iMashonkina et al.l d2007l) we estimate the NLTE 
corrections to the calcium abundances for the 6439A and 
4454A lines. For the former line, we find a minor abundance 
correction of +0.05 dex for the TO stars and no correction at 
all for the bRGB stars, while the latter line points to +0.13 dex 
at the hotter end and +0.10 dex at the cooler end. We would 
therefore expect that the trend in calcium abundance resulting 
from NLTE analysis of these lines is almost flat. Considering 
the possibility of differential NLTE abundance corrections of 
this size, it is clear that our inferred trend as such is insignificant. 



4.4. Titanium 

The titanium abundances are based on twelve Till lines and 
three Til lines, all except two (Till 5 188 A and 5 226 A) situ- 
ated in the spectral range 4 390 - 4 540 A. LTE line formation 
is believed to be a fair approximation for Ti II, since this is the 
dominant ionization stage. For neutral titanium one would on the 
other hand not expect LTE to be valid, but no thorough NLTE 
analysis exists in the literature covering this wavelength range 
for metal-poor stars. In any case, singly ionised titanium lines 
dominate neutral ones in our abundance analysis, both in terms 
of number and line strengths, and we therefore assume NLTE 
corrections to be of lesser importance. 

On the absolute scale, the Ti abundances are in good agreement 
with the results from Paper I, which are based on Till 5 188 A 
and 5 226 A. We obtain abundances ranging from 3.03 dex at 
6 250 K to 3.00 dex at 5 450 K, and the abundances from Paper I 
fall within our 1 cr limit of 0.04 dex from these values. The abun- 
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dance difference we find between TO stars and bRGB stars is 
-0.03 dex, whereas Paper I report a positive trend of 0.02 dex. 
Due to the large influence from Ti II lines, log en has a positive 
reaction to a rise in surface gravity contrary to the other elements 
and is also highly insensitive to changes in effective temperature 
at the hotter end. Lowering the surface gravity of the TO stars 
by 0.15 dex, or alternatively raising the bRGB surface gravity by 
0.21 dex, would result in a flat trend. 



4.5. Iron 

We include 51 Fel lines and 9 Fell lines in our analysis. 
Considering that the latter are comparatively weak, we estimate 
that lines from Fel bear higher weight. We find a trend of 
0.13 dex between TO and bRGB stars, again slightly larger than 
the 0.10 dex reported in Paper I. The average iron abundance of 
the sample ranges from 5.09 dex to 5.22 dex, corresponding to 
[Fe/H] = -2.41 to -2.28, using loge Fe0 = 7.5. To flatten the 
trend, one would need to raise the TO star by 3 50 - 390 K 
or low er r eff at the bRGB by 190 - 230 K. Following iKorn et al l 
(2003), Paper I find a minor NLTE abundance correction for 
neutral iron at the level of 0.03-0.05 dex in TO stars as well as in 
RGB stars and practically no differential NLTE effects between 
the groups. The small size of the corrections is likely due to the 
adoption of high rates of collisions with neutral hydrogen, and 
we do not exclude that NLTE corrections may adjust our iron 
abundances significantly upward. 

Neglecting differential NLTE influences and other possible 
modelling deficits, the significance level of the identified trend 
in iron abundance depends only on errors in relative stellar 
parameters. An estimate of the error in AlogeFe between TO 
and bRGB stars can be obtained by adding the individual 
errors listed in Table 3 in quadrature for the hottest and coolest 
stars, respectively. Assuming we have constrained the range 
of the effective-temperature scale to better than 100K, the 
surface-gravity scale within 0.1 dex, and the microturbulence 
scale within 0.1 km s _1 , the total error in abundance is 0.050 dex 
at r eff = 6250 K and 0.080 dex at T eff = 5450 K. A lower limit 
on the significance of the trend is thus 0.13/0.080 = 1.6sigma, 
corresponding to approximately 90%. 



4.6. Comparison to other studies 

In general our results for the TO stars agree well with the LTE 
abundances ob t ained for the sample included in the study by 
iThevenin et al. (2001), whilst their NLTE values are signifi- 
cantly higher. Taking the average of their LTE values (given in 
their Table 4) results in abundances that lie up to 0. 1 dex above 
the findings in this work. Similar difference is seen between our 
iron abundances and the LTE values obtained in Paper I. There 
is a larger offset, ~0.15dex, between the LTE abundances for 
magnesium and calcium. We estimate that the small differences 
in stellar parameters between the studies can account for up to 

0. 1 dex in logarithmic abundance and that the remaining part 
originates in line selection and modelling techniques. This is 
confirmed by a re-analysis of two of the UVES targets in Paper 

1, using the same automated procedure and line list as we have 
implemented in this study. However, for the discussion of abun- 
dance trends, differences in absolute abundance on this level are 
irrelevant. 

Other values reported from medium- or high-resolution spec- 
troscopic studies of NGC 6397 are [Fe/H] = -2.03 ± 



Table 3. Sensitivity of logarithmic abundances to stellar param- 
eters. 
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T c ft 


T eff +100 K 


logg + 0.1 dex 


( + 0. 1 km s 1 


Mg 


6 250 


0.053 
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Mg 


5 450 
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Ca 


5 450 
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Fig. 9. Comparison between the abundance trends from Fig. [7] 
(here displayed with black dots) and the predictions from stellar 
structure models including atomic diffusion and turbulent mix- 
ing. The horizontal, black, solid lines represent the initial abun- 
dances of the models. The three arrows in the upper left corner of 
each plot indicate the abundance change at — 6 250 K when 
raising the effective temperature by 100K (left), raising the sur- 
face gravity by 0. 1 dex (middle), and raising the microturbulence 
values by 0.1 km/s (right). 



0.02 + 0.04 deration et al.l 2001), [Fe/H] = -2.0 + 0.0 5 



(Castilho et al. 



120001). fFe /HI = -2.21 + 0.05 (lGrattonlll982l) 
-1.88 (lGratton&Ortolanilll989l). -1.99 + .01 dMinniti et al.l 
119931) . -1.82+0.04 (ICarretta & Grattonl 19971) . We note that sev- 
eral of the mentioned studies are performed exclusively on giant 
stars. 
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5. Discussion 

Abundance variations in GCs of the kind our analysis suggests 
may be explained with stellar-structure models that allow for 
radial diffusion of chemical elements. We compare our obtained 
abundances to predictions from models of Population II stars by 
Rich ard et al.l d2002ll2005l) . which account for all the physics of 
particle transport that can be modelled from first principles (see 
Paper I, Richard et al. 2005, and references therein). 
The main free parameter in these models is the introduction 
of turbulent transport of unknown physical nature, which the 
authors conclude to be best parameterised as a (decreasing) 
function of temperature and density. Additional mixing is 
needed to meet the strong observational evidence of a flat and 
thin plateau of Li in metal-poor TO and SGB stars over a wide 
range in metallicitie s, known as the 'Spite plateau' of lithium 
(Spite & Splti U982J). Omitting turbulent mixing, the plateau 
cannot be reproduced by the models, as too much Li diffuses 
into the stellar interior. In Fig. [9] we illustrate the predictions 
from three models implementing different values of a reference 
temperature To, which controls the efficiency of the mixing. 
The models are denoted by T5.80 (lowest efficiency), T6.00 
(the best fitting efficiency to the results in Paper I), and T6.09 
(highest efficiency) and represent the range in least efficient 
mixing that is compatible with the Spite plateau of lithium. The 
absolute abundances predicted by the three models are slightly 
adjusted, individually for each element, to agree with our de- 
rived abundances in the cooler half of the effective-temperature 
range (< 5800 K), where the three models coincide. Here, 
the region mixed by turbulence is fully encompassed by the 
convective envelope and the mixing efficiency is irrelevant for 
the surface abundances. The black solid horizontal lines mark 
the initial abundances of the models (accounting for the vertical 
shifts), loge Mg = 5.56, logeca = 4.41, log en = 3.07, and 
log e Fe = 5.23, or [Fe/H] = -2.27. 

The comparison shows that the predicted increase in surface 
abundances from T s s « 5 800 K to 5 400 K is overall well 
reproduced by our obtained values. We further conclude that, 
within reasonable error bars, the best-fitting model found in 
Paper I (T6.00) is also a good choice for this sample of stars. 
Generally, it is clear that the surface iron and magnesium 
abundances predicted from models including atomic diffusion 
and mixing are more compatible with the results of our analysis 
than flat trends are. 



6. Conclusions 

We have presented a homogeneous analysis of medium-high res- 
olution spectroscopic data covering a large sample of stars, lo- 
cated between the TO point and bRGB of the globular clus- 
ter NGC 6397. The obtained iron abundances show a signifi- 
cant trend with evolutionary stage. The magnesium abundances 
also indicate the presence of a trend of similar size, but the 
high sensitivity of this element to stellar parameters combined 
with the possibility of differential NLTE corrections prevent an 
unambiguous detection. The calcium and titanium abundances 
can be reconciled with a flat behaviour. The difference in iron 
abundances between the stars at T e g^ a = 5 450 K and stars at 
T e ff,Ha = 6250K amounts to 0.13dex in logep e , a trend that is 
robust to realistic errors in stellar parameters. Raising the effec- 
tive temperatures of the hottest stars by ~ 350 K would remove 
the trend in iron abundance. Considering that the total range in 
effective temperature between the considered points is 800 K, 



this would imply an error in our r e ff-scale of 43 %, which we re- 
gard as unlikely. Alternatively, the same effect may be achieved 
by lowering the effective temperatures of the coolest stars by 
~ 200 K, correspoding to a systematic error of almost 25%. 
Even that is an improbably large correction, given the good 
agreement between spectroscopy, through profile-fitting of Ha 
and the different photometric calibrations that were exploited. 
The excitation equilibrium of Fel points to cooler r e ff-values 
overall, but on the relative scale the agreement is satisfactory. 
Overall, the results of this stu dy independently support the con- 
clusions of iKorn et al.l d2007l) . indicating that atomic diffusion 
significantly affects the surface abundances of metal-poor stars 
near the main turn-off point. 

The abundance analysis is based on traditional LTE analysis us- 
ing ID, hydrostatic model atmospheres. Although it seems that 
corrections to the obtained abundances of Fe and Mg by ID, 
NLTE modelling have no large differential impact in compar- 
ison to the sizes of the found trends, we cannot rule out that 
our results are artifacts from insufficient modelling techniques. 
More sophisticated models, implementing NLTE line formation 
in 3D, dynamical model atmospheres may provide a more def- 
inite answer. Some explorations with LTE line formation using 
3D models were performed in Paper I, indicating that the results 
are robust in this respect. 

The obtained abundance trends with effective temperature are, 
for magnesium, calcium, and iron, well reproduced by stellar- 
structure models including atomic diffusion and turbulent mix- 
ing. Especially, the good relative agreement seen at the cooler 
half of the SGB, where additional turbulent transport is irrel- 
evant, indeed lends support to the notion that atomic diffusion 
shapes the surface abundances of unevolved metal-poor stars. 
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Table 4. Stellar parameters and abundances derived for the sample stars. Targets marked with * are the targets of Korn et al. (2007). 
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4.26 


3.10 
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6210 
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1.71 
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4.26 
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0.568 


6173 
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15.888 
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5440 
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0.890 
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17 40 43.750 


53 37 17.40 


15.339 


5430 


3.36 


1.60 


5.46 


4.38 


2.98 


5.25 


0.920 


5466 


5385 


0.778 


5386 


5298 


2978 


17 39 37.960 


53 40 53.10 


15.935 


5810 


3.70 


1.70 


5.24 


4.13 


2.89 


5.00 


0.755 


5964 


5772 


0.647 


5840 


5712 


3124 


17 39 38.870 


53 38 07.20 


16.136 


6190 


3.91 


1.88 


5.31 


4.33 


3.00 


5.04 


0.698 


6173 


5917 


0.590 


6077 


5917 
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5.52 


4.35 


2.96 
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0.873 


5596 
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5471 


5377 


3602 


17 39 41 630 


53 34 5(S 90 


15.981 


5840 


3.71 
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5.47 


4.23 


2.96 


5.07 
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5821 


0.628 
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16.231 


6210 


3.93 
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16.292 
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1 79 
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3 05 
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1 7 3Q 43 9gn 
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JJ JJ U7.JU 
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5490 


3.48 


1.36 


5.51 


4.37 


2.95 


5.11 


0.835 


5705 


5579 
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5549 


5450 


3892 


1 7 39 43 720 

± / J ,/ ^ J . / Z.VJ 


53 40 50 60 


15.966 


5870 


3.73 


1.50 


5.41 


4.19 


2.97 


5.06 


0.742 


6010 


5805 


0.635 


5889 
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1 / J ~ ii . / UVJ 
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J J J" T-U-JVy 


16.263 


6300 


3.99 
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5.41 


4.29 


3.04 


5.08 


0.683 


6231 


5955 
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6160 


5988 
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JJ J J. JVJ.T^v/ 


15.621 


5540 


3.52 
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5.47 
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J vy vy vy™ J 
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J J JZ. JJ.UU 
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16.013 
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4.28 
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5.05 


0.723 


6081 
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0.616 


5966 


5822 
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53 33 20 70 

J J J J — V./ . / V / 


15.514 


5440 


3.39 
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5.47 
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5.20 


0.900 
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5338 
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53 33 48 (SO 

JJ JJ T^o.vjvy 


15.906 
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3.83 


1.87 


5.38 


4.27 


2.93 


4.94 


0.768 
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5742 


0.659 
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5672 
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S01 9fS9 
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1 6 1 37 


6710 


3 93 
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5 43 


4 34 


3 06 
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6174 


5917 
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6077 
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jy io 
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1 7 40 3 1 820 

L / ^VJ J L • U Z. VJ 


S3 34 14 30 

JJ J1 1T.JU 


16.180 


6270 


3.97 


1.92 


5.29 


4.36 


3.05 


5.10 


0.691 


6199 


5934 


0.582 


6112 


5947 
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JVJ 1J07 
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3.04 
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5959 
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J J J J JZ.. / V 
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3.85 


1.76 
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4.29 
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6122 


5882 
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1 7 40 46 740 
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53 36 49 10 
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16.016 


6260 
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4.22 


2.88 


5.15 


0.721 


6086 


5858 
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5971 
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1 7 39 22 1 40 

LI ~J J Z.Z.. X^VJ 


53 38 23 40 

JJ J(J Z.J.T^V/ 


15.831 


5730 


3.65 


1.56 


5.35 


4.29 


2.95 


5.14 


0.799 


5818 


5665 


0.689 


5680 


5569 


502740 


1 7 39 10 480 
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15.607 


5430 


3.36 


1.27 


5.49 


4.38 


3.07 


5.24 
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5587 


5486 
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53 38 51 90 
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15.801 


5530 


3.51 
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4.26 
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5636 
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4.28 
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Table 5. List of lines used in the abundance analysis, with references to the oscillator strengths used. The columns Waals. , Stark. , and Rad. list 
the van der Waals, Stark, and radiative damping data (for a description of the data sources, see Barklem et al. 2005, section 3.2, paragraph 5). 
Rad. is the logarithm (base 10) of the FWHM given in rad s . Non-zero values of Stark, and negative values of Waals. represent the logarithm 
of FWHM per unit pertuber number density at 10 000 K, given in rad s -1 cirr 3 . Positive values of the van der Waals parameter correspond to the 
notation of Ans tee & O'Maral (JT995), where a packed parameter is used for the broadening cross-section (a) for collisions by neutral hydrogen 
at 10km s -1 and its velocity dependence (a). E.g. for the Mgl 5172.684A line, Waals. = 729.238 means a = 729 in atomic units and a = 0.238. 
d/lbiuc and dA Icl i give the extent of the spectral windows used for fitting of each line, blue-wards and red-wards of the central wavelength. Zero 
values in both columns indicate that the line-contribution is considered in the same window as the one listed directly above. 
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